Spin-Polarized Electron Transport at Ferromagnet/Semiconductor Schottky Contacts 



J. D. Albrecht 

Air Force Research Laboratory, Wright-Patterson AFB, Ohio 45433 



m 
O 
O 

(N 
X> 

PL, 



> 

in 

(N 
O 
m 
O 



i 

c 

o 
o 



13 



D. L. Smith 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545 

We theoretically investigate electron spin injection and spin-polarization sensitive current detection at Schot- 
tky contacts between a ferromagnetic metal and an n-type or p-type semiconductor. We use spin-dependent 
continuity equations and transport equations at the drift-diffusion level of approximation. Spin-polarized elec- 
tron current and density in the semiconductor are described for four scenarios corresponding to the injection or 
the collection of spin polarized electrons at Schottky contacts to n-type or p-type semiconductors. The transport 
properties of the interface are described by a spin-dependent interface resistance, resulting from an interfacial 
tunneling region. The spin-dependent interface resistance is crucial for achieving spin injection or spin polar- 
ization sensitivity in these configurations. We find that the depletion region resulting from Schottky barrier 
formation at a metal/semiconductor interface is detrimental to both spin injection and spin detection. However, 
the depletion region can be tailored using a doping density profile to minimize these deleterious effects. For 
example, a heavily doped region near the interface, such as a delta-doped layer, can be used to form a sharp po- 
tential profile through which electrons tunnel to reduce the effective Schottky energy barrier that determines the 
magnitude of the depletion region. The model results indicate that efficient spin-injection and spin-polarization 
detection can be achieved in properly designed structures and can serve as a guide for the structure design. 

PACS numbers: 85.75.-d, 73.50.-h, 73.40.Qv, 73.30.+y 



I. INTRODUCTION 

Semiconductor physics is in the midst of a wide rang- 
ing exploration of physical phenomena and device concepts 
that are connected with the electron spin degree of freedom. 
Much work has focused on optical generation and detection 
of spin populations. However, most spin based device con- 
cepts require an electrical means of injecting, manipulating, 
and detecting spin-polarized electron currents. Thus it is im- 
portant to understand the fundamental physics of electron 
spin transport in the main structural components that make 
up semiconductor devices. The Schottky contact is an es- 
sential semiconductor device component. Schottky contacts 
form at most metal/semiconductor interfaces. Electrical spin 
injection and detection schemes often involve ferromagnetic 
metal/semiconductor interfaces with Schottky contacts so it is 
important to understand spin dependent electron transport at 
these structures. 

Presently, theories for the injection or detection of 
spin-polarized electron currents at a metallic ferromag- 
net/nonmagnetic semiconductor interface have been formu- 
lated in the spirit of transport at a ferromagnetic metal/normal 
metal interface. The description of spin transport is incorpo- 
rated using variations on a spin diffusion equation. 1 In these 
approaches^iili^ii the semiconductor is described as a poorly 
conducting metal, in the sense that the the carrier density and 
thus the conductivity of the semiconductor is taken to be spa- 
tially uniform. Important insights gained through these mod- 
els include: (z) the large conductivity mismatch between a 
highly conductive metal and a comparatively weakly conduc- 
tive semiconductor is a major obstacle to spin injection; and 
(if) a spin selective interface resistance can be of great benefit 
to efficient spin injection. A major drawback to a spin de- 
vice physics model based on such uniform conductivity treat- 



ments is that they do not describe the underlying electronic 
properties, the currents and potentials of real semiconductor 
structures. An obvious example of a problem with uniform 
conductivity models for metal/semiconductor Schottky con- 
tacts is that they yield the symmetric, linear current-voltage 
characteristics of resistors rather than the rectifying character- 
istics of diodes. An initial study of spin injection including 
the effects of band-bending in a depletion region at an n-type 
Schottky contact showed that the depletion region can have 
an important effect on spin transport and that a device-physics 
approach to the theory of spin-contacts is necessary^ 

Experimentally, spin dependent transport has been investi- 
gated at interfaces consisting of a ferromagnetic metal 9 42ill or 
a heavily doped spin-polarized semiconductor 1213 contact and 
a non-magnetic semiconductor. Both spin injection, in which 
the electron flux flows from the spin polarized contact into the 
nonmagnetic semiconductor, and spin detection, in which the 
electron flux flows from the nonmagnetic semiconductor into 
the spin polarized contact, have been considered. In the spin 
injection measurements, detection of spin-polarized injection 
is often made using a spin-LED (light-emitting diode) con- 
figuration. In these experiments, electrons are injected into 
an n-type semiconductor from a spin polarized contact and 
are subsequently transported to a detection region, typically a 
quantum well, where they recombine with unpolarized holes 
transported from an adjacent p-type doped region. Because of 
the optical selection rules in III-V semiconductors the relative 
intensity of right- and left- circularly polarized luminescence 
gives a measure of the spin-polarization of the electron density 
in the recombination region. In spin detection measurements, 
spin polarized electrons are often optically generated in III-V 
semiconductors, and a spin dependent voltage signal is sought 
as the electron flux is transported into a spin polarized contact. 

In this paper, we theoretically investigate spin-polarized 
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electron current at ferromagnetic metal/semiconductor Schot- 
tky contacts. We systematically treat the semiconductor de- 
vice operation and the spin physics at the same level of ap- 
proximation. We consider both n-type and p-type Schottky 
contacts with current flow corresponding to either forward and 
reverse bias. We first treat the overall electrostatics of the sys- 
tem and subsequently solve charge and spin continuity equa- 
tions. We use a drift-diffusion transport model to describe the 
charge and spin currents. The drift-diffusion transport model 
is a strong scattering approximation appropriate for relatively 
high temperatures, such as room temperature. It is the ap- 
proach used to describe most semiconductor device operation. 
Here we extend this this approach to describe spin dependent 
transport at Schottky contacts. We find that the depletion re- 
gion associated with a Schottky energy barrier can have a very 
strong effect on spin-polarized electron transport at ferromag- 
netic metal/semiconductor contacts. A large Schottky barrier 
is detrimental to spin injection and can also hinder spin detec- 
tion. The model suggest structure design strategies for reduc- 
ing the detrimental effects of the Schottky energy barrier. 

The paper is organized in the following way: in Sec. II 
we describe the model, in Sec. Ill we present our numerical 
results and in Sec. IV we summarize our conclusions. Calcu- 
lational details are included in the appendices. 



II. DESCRIPTION OF MODEL 

When a metal/semiconductor interface is formed the Fermi 
energy is usually pinned within the energy gap of the semi- 
conductor. The position of the the semiconductor valence and 
conduction bands, relative to Fermi energy, at the interface 
does not depend strongly on the bulk doping of the semicon- 
ductor or on which metal is used to make the contact. For 
a given semiconductor, this energy matching position at the 
interface is largely fixed. Generally the position of the semi- 
conductor valence and conduction bands relative to the Fermi 
energy at the interface, which depends on interfacial charge 
distribution, does not coincide with the corresponding energy 
position of the bands in the bulk of the semiconductor, which 
depends on the bulk doping level. There is a band bending re- 
gion near the interface which at zero applied bias: is depleted 
of carriers, is charged because of the background doping, and 
has a large spatially varying electric field. The Schottky en- 
ergy barrier between the pinned Fermi level and the semicon- 
ductor conduction band at the interface results in the charged 
depletion region and has important consequences on charge 
current flow at metal/semiconductor interfaces. For example, 
it leads to diode type current-voltage characteristics. Thus, it 
is not particularly surprising that this energy barrier and deple- 
tion region also have important consequences on spin current 
flow at these interfaces. 

The design of the interface is central to spin injection and 
detection structures. In particular, spin dependent interface 
resistance resulting from spin dependent tunnel barriers have 
been argued to be essential for effective spin injection or de- 
tection at metal semiconductor interfaces. 4 - 5 - 6 Possible spin- 
selective interface resistance layers, formed from thin mag- 



netic insulators, have been experimentally investigated by 
Motsnyi, et al»A2 In other work, Hanbicki, et al., have investi- 
gated Schottky barriers with heavy doping near the interface 
to study structures in which current is dominated by tunnel- 
ing for spin injection. 15 These results are promising for the 
realization of future electron spin based device designs. In- 
terfacial spin-flip scattering, which would be detrimental to 
spin injection or detection structures, is possible. 16 Structures 
should be designed to minimize this process. 

We consider four scenarios corresponding to the injection 
or collection of spin polarized electron current at Schottky 
contacts to n-type or p-type semiconductors. The four cases 
are schematically shown in Fig.^ Panel (a) of Fig.^iUus- 
trates the case of spin injection into a n-type semiconductor. 
The diode formed by the Schottky contact is in reverse bias 
and the electron flux is from the ferromagnetic metal on the 
left into the semiconductor on the right. A heavily doped re- 
gion near the interface, as illustrated by the doping profile in 
the lower part of panel (a), can be designed to form a sharp 
potential profile through which electrons tunnel. The heav- 
ily doped region reduces the effective Schottky energy barrier 
that determines the properties of the depletion region^ The 
total barrier e</>{, is divided into two parts, a tunneling region 
with barrier height e0 f and an effective Schottky barrier height 
eVbi. The potential drop in the depletion region consists of the 
effective Schottky barrier height plus the applied reverse bias 
eVfl. Two parameters of the tunneling region, its tunneling 
resistance and the magnitude of the reduction of the effective 
Schottky barrier, can be separately controlled by the parame- 
ters of the doping profile, for example the height and width of 
the heavily doped region. 

Panel (b) of Fig. ^ illustrates the case of spin-polarization 
sensitive current detection at a Schottky contact between a 
ferromagnetic metal and an n-type semiconductor. A spin- 
polarized electron flux is incident from the semiconductor and 
the Schottky diode is in forward bias. In a typical experimen- 
tal situation, the structure is held under constant current bias 
and a change in voltage signal is sought when the the polarity 
of the spin-polarized incident current is reversed. There may 
be a heavily doped region near the interface as in panel (a). 

Panel (c) of Fig.^iUustrates the case of electron spin injec- 
tion into a p-type semiconductor. The p-type Schottky diode 
is in strong forward bias. There is an insulating tunneling bar- 
rier at the interface that limits the hole current, which nonethe- 
less can be considerable. There is a hole accumulation region 
in the semiconductor near the interface. The minority carrier 
electron flux is from the ferromagnetic metal into the semi- 
conductor. This structure can be interesting for characterizing 
the spin dependent transport properties of the tunneling bar- 
rier. 

Panel (d) of Fig. ^ illustrates the case of spin-polarization 
sensitive current detection at a Schottky contact between a fer- 
romagnetic metal and a p-type semiconductor. The p-type 
Schottky diode is under zero or small (either forward or re- 
verse) bias. Spin polarized electrons are optically generated 
by absorption of circularly polarized light. In a typical ex- 
perimental situation, the structure is held under constant cur- 
rent bias and a fixed incident optical intensity and a change 




in voltage signal is sought when the the polarity of the circu- 
larly polarized incident light is reversed. This is the essen- 
tially the same structure as in panel (c), except under differ- 
ent bias and optical excitation conditions. These two experi- 
mental configurations can be used together to characterize the 
spin transport properties of a tunneling barrier at the ferro- 



magnet/semiconductor interface. 

We describe the ferromagnetic metal and the interface us- 
ing a spin dependent drift-diffusion equation, a spin diffu- 
sion equation, and spin dependent interface conductances as 
in Ref. 4, The drift-diffusion equation describing current flow 
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in the ferromagnetic metal is 
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Here j v is the current density due to electrons of spin type 
?7(=T,i), c j; is the conductivity for electrons of that spin type, 
H v is the corresponding electro-chemical potential, e is the 
magnitude of the electron charge and x is position. Eq. Q 
assumes rapid wave vector randomizing scattering events, so 
that electrons of the same spin stay in local quasi-thermal 
equilibrium with each other. However, spin-flip scattering can 
be comparatively slow so that electrons of different spin may 
be driven out of local quasi-thermal equilibrium by, for exam- 
ple, an applied current density. When electrons with different 
spins are driven out of local quasi-thermal equilibrium, so that 
/if is not equal to at some point in space, spin relaxation 
away from that spatial point is described by a diffusion equa- 
tion 
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Here A c is the spin diffusion length in the metallic con- 
tact and we use the notation /i±=^j±/ij. At the con- 
tact/semiconductor interface, electrons of different spin can be 
driven out of quasi-thermal equilibrium by current flow. Far 
from the interface, as ir^±oo, the electrochemical potential 
difference vanishes fj,-— >0. The total steady state current den- 
sity is a constant function of position. We assume no strong 
spin flip scattering at the interface so that the individual cur- 
rent components for the two spin types are continuous at the 
interface. Current flow at the interface is described using an 
interface resistance 



3n = 



eR r , 



(3) 



where j° is the current density at the interface, R v is the in- 
terface resistance, and A/i^ is an interfacial discontinuity in 
electro-chemical potential for electrons of spin type r/. If the 
interface resistance is zero, the electro-chemical potentials are 
continuous at the interface whereas for nonzero values of R v a 
discontinuity in /i v can develop at the interface. For notational 
ease, we set the variables j±=j^ji- In the contact where the 
hole current is zero, the total current j=j+- We take the elec- 
tron density as a function of position to be fixed, independent 
of the current density j, in the contact. The total conductiv- 
ity of the contact is then independent of position and current 
density. It is convenient to define a contact polarization vari- 
able a c by (T^—a c a c or <tj_ = (1— a c )a c where a c is the total 
contact conductivity. 

We take the contact on the left (x<0) and the semiconduc- 
tor on the right (a;>0) of the interface located at x = 0, as in 
Fig-E so that the current density is negative for electron in- 
jection into the semiconductor. Solving Eq. (0 with the stated 
boundary conditions gives 
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Quantities evaluated at the interface approached from the con- 
tact and semiconductor are indicated by the superscripts 0~ 
and + , respectively. From Eqs. and (0} we find 
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The total current in the semiconductor is 
i = JT+il+.?p = J++ip where j p is the hole current den- 
sity. For n-type Schottky barriers j p =Q, but not necessarily 
for the p-type structures. It is convenient to define (3 as the 
fraction of the electron current carried by spin up electrons 
(3=j]/j+- We assume that there is not strong spin-flip 
(2) scattering at the interface so that j . =j^_ . The total current 
j is continuous at the interface. The interface resistance 
conditions, Eq. (0 then lead to the interface matching 
conditions, 
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(7) 
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for x < 0. 



(4) 



These matching conditions apply for each of the four cases 
illustrated in Fig. Q. 

The semiconductor near the interface is either depleted, as 
shown in panels (a), (b) and (d) in Fig. [2 or accumulated as 
shown in panel (c) in Fig.^ We input the drop in electrostatic 
potential between the semiconductor side of the interface and 
the edge of the depletion or accumulation region. A tunnel- 
ing region, such as is illustrated in panel (a) in Fig. is de- 
scribed by the interface resistances i?| and R±, and is taken 
to have negligible width. The semiconductor side of the in- 
terface starts at the right of the tunneling region. From the 
input potential drop we calculate: the current density; the net 
bias voltage, which may include a contribution from the inter- 
face resistance; and the electrostatic profile. For the depleted 
cases we use the usual depletion approximation to describe 
the electrostatics in the semiconductor. For the accumulated 
case, we assume that hole current is limited by an interfacial 
barrier, take a constant hole quasi-Fermi energy in the semi- 
conductor as shown in Fig. [He), an d solve Poisson's equa- 
tion self-consistently to determine the electrostatic potential 
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in the accumulation region. Details of the electrostatics are 
described in Appendix A. 

In the semiconductor, electron and hole currents satisfy 
continuity equations, 



^L = -e(g v -r v ), 



and 



dip 
dx 



(10) 
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where g is a generation rate and r a recombination rate. For 
spin polarized electrons there is a contribution to the recom- 
bination rate from both spin flip scattering and electron-hole 
recombination, 



/ r 



(12) 



where r r and r s are the recombination and spin flip times, 
respectively, and an analogous expression applies for r± . We 
use a drift-diffusion approximation to describe electron and 
hole currents, 
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where p n M is the electron (hole) mobility, rii is the intrinsic 
carrier concentration and <p is the electrostatic potential. The 
carrier densities are given by 
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and 



P = n t e 



-(e<P+fj, p )/kT 



(15) 
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It is convenient to go into a representation describing the 
electron charge and spin degrees of freedom and we define 



so that 



and 
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dx 
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The corresponding generation and recombination rates are 

g±=g-f±gi, r + =n + /r r , and r_=ra_(2/T s +l/7y). The con- 
tinuity equations become 



dj± 
dx 



-e(g± - r±) . 



(20) 



Substituting the drift-diffusion form into the continuity equa- 
tion gives a transport equation for fl±, 
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and 



where g±=g±/\fl n (kT/e){n i /2)], A 2 + ={kT/e)fl n . 
A 2 _ = (kT/e)fi n (2/T s + l/r r ) _1 . An analogous equation 
holds for holes. Analytic solutions for f2± are discussed in 
Appendix B. 

Boundary conditions at the semiconductor side of the inter- 
face, x=Q + , and at the depletion edge, x—w, are used to de- 
termine the two matching coefficients (see Appendix B) that 
appear in the solutions of Eq. fl21i . Details of the boundary 
conditions differ somewhat for the individual cases and will 
be specified in the discussion of these cases. For the charge 
degree of freedom, we use interface recombination boundary 
conditions at the semiconductor side of the interface, 



if 



sr (»f-n£) 



(22) 



where v sr is the surface recombination velocity and n° eq is 
the equilibrium electron density at the semiconductor side of 
the interface. For the n-type semiconductor cases we set the 
electron density at the depletion edge (x=w) equal to the bulk 
doping density so that that the material becomes charge neu- 
tral at this point n + (w)=Nd where N4 is the bulk doping den- 
sity. From the definition of f2±, we see that 



/i_ = 2fcTtanir 1 ( — 

i l-L- 



(23) 



Combined with Eq. this gives a boundary condition for fi- 
at the semiconductor side of the interface. It is often useful to 
write 



0_ 



dlnQ 4 

dx 
dlnQ- 

dx 



(2/3-1) 



and 



(2/3-1) = 



J+ 



dx 
dx 



(24) 



(25) 



This form can be useful because (2/3—1) can become the un- 
known in the matching condition of Eq. (|8j. In the doped 
material beyond the depletion region the electric field is small 
and spatially uniform. In the usual treatment of current flow 
in Schottky diodes this small field is neglected. For most of 
the cases there is no generation term in Eq. ( 12 II and we only 
need the homogenous solution for 57 _ 



ft-(x > w) = Cl_(w) exp[(w — x)/£-i 



where 
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Here, E is the uniform electric field in the doped material be- 
yond the depletion region. Substituting into Eq. dl9> we ob- 
tain the more familiar notation of the "drift-diffusion" frame- 
work. In the bulk, rt_ relaxes according to 



1.0 



ri-(x > w) — ri—(w) exp[(ui — x)jlJ\ 



(28) 



where l_ reflects the field modification of the diffusion length 
at constant carrier density. 2-7 Matching the continuity of f2_ 
and its spatial derivative, the current from Eq. Q, at the edge 
of the depletion region give the final boundary conditions. 



III. CALCULATED RESULTS OF MODEL CONTACTS 

We discuss results for each of the four cases shown in Fig. 
[^sequentially. 

A. Injection at an n-type contact. 

This section discusses results for electron spin injection at 
an n-type Schottky contact as shown in Fig.^a). The effec- 
tive Schottky barrier of the of contact can be varied by chang- 
ing the doping profile in the semiconductor near the interface 
as shown at the bottom of the panel. We consider a heavily 
doped region near the interface that creates a narrow tunnel- 
ing region. The effect of the narrow doping region is to re- 
duce the effective Schottky barrier energy with the associated 
reduction in the depletion width. This approach to tailoring 
effective Schottky barriers is well established in semiconduc- 
tor device applications^ 

For this case, there is no hole current or optical generation. 
Consequently, we solve the homogeneous form of Eq. (12 1 1 
matched to boundary conditions at the interface given by Eqs. 
(I7I9> and at the depletion edge by Eqs. (I26> and Q. The tun- 
neling regions are parameterized using interface resistances 
and reduced effective Schottky barriers. The contact is metal- 
lic with resistance equal to 10 _5 f2 cm, polarization a c =0.9 
(80% polarized) and spin diffusion length equal to A c =100 
nm. The n-type semiconductor has an electron mobility of 
/x„=5000 cm 2 /(V s) and a spin diffusion length equal to 1.0 
//m. The diode characteristic from Eq. 1221 is determined us- 
ing w sr =10 7 cm/s. 

In Figs. E]|3] andH we show calculations of spin injection 
through a depleted n-type Schottky contact at T=300K. In Fig. 
13 we show the effect of the effective Schottky barrier on spin 
injection. In panel (a), the current spin polarization as a func- 
tion of position is plotted for three effective Schottky energy 
barriers, as labelled in the figure. For each barrier energy, the 
structure is biased to operate at 90% of the reverse-saturation 
current. The calculation shows clearly that the presence of an 
energy barrier degrades the performance of the spin injecting 
structure, and that the dependence on barrier height is strong. 
In panel (b), the corresponding electrochemical potential dif- 
ferences are plotted for each structure. We see from this panel 
that the origin of the splitting in electrochemical potentials (di- 
rectly related to polarization) is from the interface resistance. 
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FIG. 2: Calculated spin injection properties for an n-type Schottky 
contact as shown in Fig. 0a). (a) Current spin polarization (J- / 
and (b) electrochemical potential difference for spin-up and spin- 
down electrons as a function of position for various values of the 
effective Schottky barrier. The edge of the depletion region is indi- 
cated by x on the curves of part (a). 



If the interface resistance is lowered or if the i?| / R [ ratio ap- 
proaches unity, then the injection properties of the structure 
degrade. Some specifics regarding the barrier lowering and 
interface resistance for the n-type injector are discussed Ref. 

s 

The electron spin density polarization can be examined in 
the presence of the electron density profile. For the same con- 
ditions used in Figs.|2j the spin polarization of the local elec- 
tron density is shown in panel (a) of Fig. [3] The total density 
is shown in panel (b) of Fig. [3] for comparison. The polar- 
ized current may persist deeper into the semiconductor than 
its ability to spin polarize the local electron gas. A spin polar- 
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FIG. 3: Effect on the local spin populations of injection at an n-type 
contact as shown in Fig. 0a). (a) Electron density spin polarization 
(n_/n+); and (b) total electron density (n+) as a function of posi- 
tion for various values of the effective Schottky barrier. The edge of 
the depletion region is indicated by the x on the curves of part (a). 



ized current may be established in the semiconductor without 
strongly perturbing the spin polarization of a background of 
free carriers. However, for optical detection, such as the spin- 
LED, the signal is proportional to the spin polarization of the 
local density and not of the current. 

In the calculations presented in Fig.|2]and Fig.[3]the effect 
of an electric field in the doped region outside the depletion 
region was neglected. This is the usual approximation in de- 
scribing the electrical properties of Schottky diodes. It is rea- 
sonable because the doped region outside the depletion region 
is conductive and the current flow is limited by the depletion 
region. In Fig. [4] we compare a calculation of the current po- 
larization as a function of position neglecting the electric field 




N D =5x 10 16 cm" 3 
R i =2R T =10" 3 ficm 2 



50 100 
position (nm) 



FIG. 4: Effect of an electric field beyond the depletion region. The 
solid curve shows current polarization as a function of position in- 
cluding a greatly exaggerated electric field in the doped region be- 
yond the depletion region whereas the dashed curve assumes no elec- 
tric field outside the depletion region (as in Fig. |2j for Vbi=0.1V. 
The field used for the solid line was chosen to give £- = 10/j,m, or 
70 times the value determined from drift in the bulk. If the field 
determined by en+(T n and j is used, the result is essentially indistin- 
guishable from that using zero field. 



in the doped region with one which includes a greatly exag- 
gerated value for the electric field outside the depletion region 
for an effective barrier height of 0. 1 eV. The field used for the 
solid line was chosen to give £_ = 10/im and is 70 times that 
determined by the conductance of this region and the injected 
current density. If a field determined by the conductance and 
current density (E—j/en+p) is used, the result is essentially 
indistinguishable from that using zero field. (For the calcu- 
lation in Fig. [4] that field is —37 V/cm.) The figure shows 
that for a Schottky structure with a significant effective bar- 
rier height the electric field in the doped region outside the 
depletion region has little effect on the spin injection proper- 
ties of the structure. The reason for this is that the matching 
conditions on the currents and electrochemical potentials are 
at the interface between the metal and the depleted region of 
the semiconductor where the concentration of electrons is ex- 
ceedingly small. This is to be contrasted with uniform con- 
ductivity models where the electron concentration is the same 
up to the interface so that d/j, v /dx is driven by the electric 
field on the semiconductor side^ii 

The depleted region that occurs at at a Schottky contact 
is seen to be detrimental to spin injection at a ferromagnetic 
metal/ semiconductor interface. The problem arises because 
injection is into a very low resistance region of the semicon- 
ductor that is depleted of carriers. However, the depletion re- 
gion can be tailored using a doping density profile to minimize 
these deleterious effects. For example, a heavily doped region 
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FIG. 5: Current spin polarization in the semiconductor for an n-type Schottky structure operating in detection mode as shown in Fig. 0b). 
Results are plotted for various A s at fixed Vy. The electron flux, incident right-to-left at the depletion edge, corresponds to a current equal to 
5.0 A cm -2 , (a) Incident polarization j_ (w) /ji+=0.6 and zero interface resistance, (b) j-(w)/j+=— 0.6 and zero interface resistance, (c) 
j- (w) /j+=0.6 with fixed interface resistance, (d) j- (w) /j+ =— 0.6 with fixed interface resistance. 



near the interface, such as a delta-doped layer, can be used to 
form a sharp potential profile through which electrons tunnel 
to reduce the effective Schottky energy barrier that determines 
the magnitude of the depletion region. 

B. Detection at an n-type contact. 

This section describes results for spin detection at n-type 
Schottky structures as shown in Fig.^b). The spin detection 
case is similar to that of spin injection discussed in the pre- 
vious section but with a modification of the boundary condi- 
tions. The boundary condition on the current spin polarization 
at x=w is an incident polarized current j_ (w). We consider a 
constant total current density of 5.0 A cm~ 2 and seek a volt- 



age signal as the polarity of the spin polarization is reversed. 
The contact is metallic and has the same properties as in the 
previous section. 

The calculated current polarization as a function of position 
within the depletion region for detector operation is shown in 
Figs. [5] and |6] Results are calculated for 60% incident spin 
polarization. In Fig. [5] we examine the dependence of the 
spin polarized current on A s for fixed effective Schottky bar- 
rier. Since the effective Schottky barrier is fixed the depletion 
width is the same for the various cases. Long spin relaxation 
times result in larger spin polarizations at the interface. In 
panels (a) and (b) of Fig.[3]the current polarization behavior is 
almost identical for up and down incident currents when the 
interface resistance is zero. This shows that the presence of a 
polarized contact material (polarized spin up) has little impact 
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FIG. 6: Current spin polarization in the semiconductor for an n-type Schottky structure operating in detection mode as shown in Fig. 0b). 
Results are plotted for various Vbi at fixed A s . The electron flux, incident right-to-left at the depletion edge, corresponds to a current equal to 
5.0 A cm -2 , (a) Incident polarization j_ (w) /j+=0.6 and zero interface resistance, (b) j- [w)/j+=— 0.6 and zero interface resistance, (c) 
j- (w)/j+— 0.6 with fixed interface resistance, (d) j- (w) /j+=— 0.6 with fixed interface resistance. 



on the spin polarized current in the semiconductor. In panels 
(c) and (d) the currents are calculated including interface re- 
sistance and show a strong asymmetry owing to the mismatch 
in i?| and R± . 

In Fig. [5] we examine the dependence of the current polar- 
ization in the depletion region on Vu at fixed A s . The bias 
conditions have been adjusted to give the same total current 
density for all cases. Structures with different barrier heights 
have different depletion widths. A main point seen from Fig. 
|6]is that large effective barrier energies result in small spin po- 
larizations at the interface. In panels (a) and (b) of Fig.[6]the 
current polarization behavior is very close for up and down in- 
cident currents when the interface resistance is zero. In panels 
(c) and (d) the currents are calculated including interface re- 
sistance and show a strong asymmetry owing to the mismatch 



in i? j and i? j . 

We consider an n-type Schottky detector structure at a con- 
stant total current density. To fix the total current in the struc- 
ture, the forward bias (Vf in Fig.^b)) is tuned for each struc- 
ture. The detected signal is the change in voltage at fixed cur- 
rent density when the spin polarization of the incident current 
is reversed. This voltage difference is obtained by integrating 
the electrochemical potential over position from — oo to +oo 
and taking the difference for two incident spin polarizations. 
After cancellations, the surviving terms yield a voltage differ- 
ence of 

AV = A ( x i± k -R l+ Ac (1 ~ 2ac \ 1 (29) 
\j+ J 4 [ ' a c a c (l-a c )\ 
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FIG. 7: Voltage signal as a function of effective Schottky barrier 
energy for j_/j+=±0.6 and incident current equal to 5.0 A cm -2 . 
Results are shown for three values of the spin diffusion length. 



where A indicates the difference between the quantities for 
opposite signs of the incident spin polarization at x=w. 

^From Eq. j29t we see that without a spin dependent in- 
terface resistance (R^=Ri—>Q) and for a highly conductive 
contact (er c ^oo) no significant voltage difference can be es- 
tablished. Using metallic contacts (a c —10 5 13 cm -1 ) with 
contact spin diffusion lengths less than a micron (A c =10 -5 
cm) and currents (~1 A cm~ 2 ) corresponding to low biasing 
conditions the calculated detected voltage differences are neg- 
ligibly small (<10 -10 V). We conclude that Schottky contacts 
without a spin-selective tunnel barrier will not be useful as 
spin polarized current detectors and therefore we concentrate 
on Schottky structures containing a spin selective interface re- 
sistance. In Fig.0we show calculated voltage differences as 
a function of effective Schottky barrier for three values of the 
spin-diffusion length in the semiconductor. The detection sig- 
nal saturates for both large and small Vu values. For small 
barriers the depletion region vanishes and the incident polar- 
ized current reaches the interface and only small applied bias 
is required to establish the constant current. For large deple- 
tion widths, larger applied biases are required to keep the cur- 
rent density fixed and the resulting resulting interface current 
polarization saturates consistent with the behavior shown in 
Fig. He) and^d). 

The depleted region that occurs at at a Schottky contact is 
seen to be detrimental to spin detection at a ferromagnetic 
metal/ semiconductor interface. The problem arises because 
in these forward biased structures, electron current is driven 
by diffusion against a strong and rapidly varying electric field 
in the depletion region. As a result the effective drift-diffusion 
lengths in the depletion region can become rather short lead- 
ing to strong spin relaxation. As for the electron injection 



structures the depletion region can be tailored using a doping 
density profile to minimize these deleterious effects. 



C. Injection at a p-type contact. 

This section describes results for spin polarized electron in- 
jection from a ferromagnetic contact into a forward biased p- 
type Schottky diode. The hole current is limited by a barrier 
at the interface but is most likely larger that the minority elec- 
tron injection current. The output signal is the ratio of right to 
left circularly polarized light emitted from the semiconductor 
when the injected electrons recombine radiatively with unpo- 
larized holes in the p-type material. 

The electrostatic treatment for the p-type structure is de- 
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accumulated p-type Schottky diode in strong forward bias as in Fig. 
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FIG. 9: Integrated emitted photon polarization 

— CT_)/(cr++cr_) as a function of electron spin-flip scat- 
tering rate coefficient for fixed injected current density from a 
ferromagnetic contact (a c =0.9) into the accumulated p-type contact 
shown in Fig. [5] Results are shown for injection efficiencies of 
j+/j =0.5 (solid curves) and 0.1 (dashed curves). 



left circularly polarized light (a±) is proportional to the local 
electron density and is given by 



r 

Jo 



.dx 



J °° n+dx 



A 2 _ 
Ai 



(2/3 0+ 



(30) 



The spin polarization (2/3 0+ — 1) of the electron current at the 
interface is obtained from the numerical solution of Eq. 
The integration is straightforward because the electron con- 
centrations are related to the current through the continuity 
equations, Eqs. dl0ll2l >. 

The calculated optical polarizations for injection into p- 
type material are shown in Fig.|9]as a function of the spin-flip 
scattering rate coefficient. The curves are calculated for the 
same electron current density at x=0. We consider insulat- 
ing barriers that result in minority carrier injection efficiencies 
equal to 50% and 10%. If used as a characterization tool, the 
structure should to be sensitive to differences in the optical po- 
larization signals in order to determine the interface resistance 
values. This is the case for the case of slow spin-flip scattering 
compared radiative recombination. If the spin-flip scattering 
rate is faster than the recombination rate in the semiconductor 
characterization by this method is will to be difficult. 



D. Detection at a p-type contact. 



scribed in Appendix A. Calculated energy band and hole con- 
centration profiles are shown in Fig. [8] and correspond to the 
semiconductor portion of Fig. ^c). An interfacial barrier is an 
important feature of the structure to prevent runaway hole cur- 
rent at the interface. As a consequence of the barrier there is 
strong accumulation of holes near the interface. The large and 
rapidly varying hole concentration requires a different treat- 
ment of the electron spin relaxation process 18 than used in the 
depletion cases considered above. To account for the presence 
of the accumulated holes, we use the local hole concentration 
to vary the recombination and spin-flip scattering rates as a 
function of position. 

We take the electron spin relaxation times to be linear 
with the local hole density and define rate coefficients R s t 
and R rec for spin-flip scattering and recombination so that 
T~ 1 =R s fXp(x) and T J T 1 —R rec xp(x). We present calcula- 
tions for a range of R s f. We take R rec —7x 10~ 10 cm 3 s _1 , a 
typical value for p-type GaAs at T=300K. 19 The hole mobility 
is/ip=500cm 2 /(Vs). 

Unlike for the depleted structures there is no closed form 
analytic solution for the spin-dependent transport equations 
for the accumulated structure. Given a numerical solution for 
p(x) from the electrostatic calculation, the solution of Eq. J21t 
is obtained by numerical integration and the shooting criterion 
that solutions be non-diverging as x-^oo. There is no optical 
generation of carriers (g±=0) so that only the homogeneous 
solution is required. 

The detected quantity in this case is the degree of circular 
polarization of the emitted light. Assuming good radiative re- 
combination efficiency, the optical polarization of right and 



The p-type detector structure shown in Fig. ^d) involves 
optical generation of carriers. The p-type Schottky diode is 
at zero or small bias (either forward or reverse). As in the 
depleted n-type detector, we compute the voltage difference 
established as the circular polarization of the light is changed 
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FIG. 10: Voltage difference for collection of spin polarized electrons 
at a depleted p-type contact. The inset shows the electron and hole 
current contributions to the current-voltage characteristics. 
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from right to left. We assume that the optical generation is 
nearly uniform over the depletion depth (i.e., that the recipro- 
cal of the absorption coefficient is small compared to the de- 
pletion width) and that the optical generation of spin-polarized 
carriers is a 3:1 or 1:3 ratio of spin-up to spin-down electrons 
according to the optical selection rules. Beyond the depletion 
edge, for x>w, the optical generation falls off as exp(— \x) 
with x=10 4 cm. Under these conditions, the transport equa- 
tion (12 II is solved analytically as described in Appendix B. 

A calculated current/voltage characteristic for the various 
current components calculated for this structure are shown in 
the inset of Fig. The properties of the p-type Schottky 
are those used in the accumulation case of the previous sec- 
tion. The barrier height is 0.7eV and the incident light power 
is set to 1 W cm~ 2 or g + «4.4xl0 22 s" 1 cm~ 3 for photons 
at the bandgap of GaAs. In this structure, the electron cur- 
rent stays nearly constant over a broad range of applied bias 
around Va=Q so we report results for zero bias. The value 
for AV, as for the n-type detector structure, is given by Eq. 
d29t except that the total electron current at the interface j + 
must be accounted for separately from the total current which 
includes a contribution from holes. This separation is straight- 
forward once solutions for £l± p have been calculated. In Fig. 
|lO|we plot the calculated voltage differences as a fucntion 
of A s for optical excitation. If the semiconductor spin life- 
time has been determined by other experimental means, this 
represents a second method of characterizing the interface re- 
sistance. Notice that even when the spin lifetime becomes 
maximally long (A s =A r ) there is a strong dependence of the 
measured voltage on the interfacial conditions. 



IV. SUMMARY 

We have presented a theoretical description of the injec- 
tion and detection of spin polarized electrons at n-type and p- 
type Schottky contacts. The presence of the depletion region 
that occurs at Schottky contacts has been shown to be detri- 
mental to both spin injection and detection in n-type Schottky 
structures. The problem for the reverse biased n-type spin in- 
jection structures arises because the injection is into a very 
low resistance region of the semiconductor that is depleted of 
carriers. The problem for the forward biased n-type spin de- 
tection structures arises because electron current is driven by 
diffusion against a strong and rapidly varying electric field in 
the depletion region. As a result the effective drift-diffusion 
lengths in the depletion region can become rather short lead- 
ing to strong spin relaxation in the depletion region. 

For both n-type injection and detection structures, the de- 
pletion region can be tailored using a doping density profile 
to minimize these deleterious effects. A heavily doped region 
near the interface, such as a delta-doped layer, can be used 
to form a sharp potential profile and this tunneling region ef- 
fectively reduces the Schottky energy barrier that determines 
the magnitude of the depletion region. The model results in- 
dicate that efficient spin-injection and spin-polarization detec- 
tion can be achieved in these n-type structures if they are prop- 
erly designed so that the effective Schottky barrier is reduced 



to less than about 0.2 eV. 

We discussed two experimental cases in which ferromag- 
netic Schottky contacts to p-type semiconductors could be 
used to characterize the spin dependent transport properties 
of interface tunnel barriers: optical detection of spin cur- 
rents injected into a strongly forward biased accumulated p- 
type semiconductor and electrical measurement of optically 
excited spin populations at zero biased p-type contacts. The 
same structure can be used under different bias and excitation 
conditions for the two experiments. 

A set of spin-labelled transport equations has been devel- 
oped, in a systematic way, that is suitable for device models 
at a basic level. We have demonstrated that the electrostatic 
and current conditions that are present in actual devices can 
lead to important consequences for spin dependent transport 
in the structures must be taken into account. With these build- 
ing blocks in place, there are clear extensions of the model 
to more complex device structures which will be the focus of 
future work. 
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APPENDIX A: SCHOTTKY CONTACT ELECTROSTATICS 

The band-bending that occurs near a Schottky contact is 
central to the discussion of spin transport in these structures. 
This appendix describes the electrostatic inputs to the solu- 
tions of the spin transport equations. 

Three of the cases, shown in Fig. [Ha), (b), and (d), involve 
a contact under small bias. In these cases, we use the depletion 
approximation where the electrostatic potential and depletion 
width are given by 

(f){x) = 0(0) ± —tot -^-x 2 0<x<w (Al) 

£ s^o 2 £ s£o 

and 



w= ^^^(V bl + V A ) (A2) 

here e is the free space permittivity, and e s is the relative 
static dielectric constant of the semiconductor. The applied 
voltage Va can be either sign but must be small enough in 
forward bias so as not to invert the semiconductor from deple- 
tion to accumulation. For an n-type contact, N is the donor 
concentration and the upper sign applies. The acceptor con- 
centration and the lower sign are used for p-type contacts. 

The inverted p-type contact shown in Fig. ^c), requires 
a numerical treatment of the electrostatics. We consider a 
case of strong forward bias in which a barrier layer limits 
hole transport at the metal/semiconductor interface and we 
treat the accumulation region in the semiconductor as a quasi- 
equilibrium system characterized by a Fermi energy for holes 
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with a valence band density of states Ny- The electrostatics 
in the accumulation region is determined by the hole density 
which is given by statistics: 



which can be simplified by substituting the analytic expres- 
sion for the interior integral. After substitution and simplifi- 
cation, we arrive at the following result: 



p(x) = N v -=^ [\ F {x)] 

where the Fermi one-half integral is given by 

A^ 



o i + e x -y 



dx 



(A3) 



(A4) 



and 



. , , s- E v (x) , . Ep — Ey (x) 
X (x) = Xf (x) = — w . (A5) 



kT 



kT 



Far from the interface the hole density is equal to the acceptor 
doping level p(oo)—N A so the local hole concentration can 
be written as 



P (x) = N A 



^ [Mg] 
T k [A F (oc)] ' 



(A6) 



To solve Poisson's equation, we compute the charge density 
by subtracting the background density of ionized accepters 
which yields 

p{x) = e\p (x) - N A ] 

Tx \\p{x)\ — T\ \Xp(oo)] 

= eN A ^ * . (A7) 

Ti [A F (oo)] 

These algebraic steps are made so that the numerator can be 
rearranged to give an integral with a closed-form solution. 
First we rearrange the difference of Fermi integrals in the nu- 
merator as: 



r h [\p{x)}-T.[Xp{oo)] 



A 1 / 2 



where 



A 



J _ e \F(oo)-\F{x) 

e \p(oo) — \ _|_ gAir(oo) — Af(i) 
Ey(.)-E r [oo) 

1 — e 

e A F (oo)-A _|_ g— jpjr-* 



AdX (A8) 



(A9) 



We substitute the charge density into the integral form of Pois- 
son's equation for the electric field F. This process begins 
with Poisson's equation in the form 

F(oo)=0 rEv(oo) / p \ 

FdF= / ti^LldEv (A10) 

JE v (x) £ s£o 



>F(x) 

and results in 

F 2 (z) 



e s e T]_ [Af(oo)] 

\l/2 / fE v (oo) 



la 1 + exp(A- Aoo) W Bv(:c) 



AiiV dA (All) 



F 2 (x) 



2kTN A 



e s e Ti_ [Af(oo)] 

A l/2 



1 



qg + (1 + 9) In 



1 + qe~ 
1 + 



dA (A12) 



where q and (7 are defined by 



Ey (00) — Ey (x) , , 

g = fcy a = exp(A- Aoo). (A13) 

An electrostatic profile is computed by first obtaining the 
electric field at x=0. Given the energy barrier (fixes q at x=Q) 
and doping density (fixes the fermi level relative to the band 
edge in the bulk), we can obtain F(Q) with one numerical 
integration of Eq. JA12I >. The valence band, electric field, and 
density profiles are then generated by integrating Eqs. jAlCH 
and JA12t forward in x and evaluating Eq. JA7> at each spatial 
step. A sample result of this process is plotted in Fig. [8] 



APPENDIX B: SOLUTIONS TO THE TRANSPORT 
EQUATION IN THE DEPLETION REGION 

In this appendix we present analytic solutions for the trans- 
port equation (12 1 1 when the depletion approximation is used 
for the electrostatics of the structure. 

In the depletion mode cases, we have a second order differ- 
ential equation in the depletion region that is solved subject to 
boundary conditions at the metal/semiconductor interface and 
at the depletion region edge in the semiconductor. A typical 
equation, including the possibility of optical generation, can 
be written as 



d 2 n 

dx 2 



e d(f> dQ 
kT dx dx 



1 



A 2 



(Bl) 



Within the depletion approximation, we have a quadratic form 
of the electrostatic potential 



e4> 1 2 e(p 

— = —ax + ox H 

kT 2 kT 



The particular solution fl p is 



G 



P a+ 1 



■e kT 



(B2) 



(B3) 



A 2 



and the homogeneous differential equation can be written as 

(B4) 



d 2 n . L .dn 1 _ n 

+ (ax + b) — rfl = 0. 

dx 2 dx A 2 



Solutions results from a change of variables that transform 
the homogeneous differential equation to the confluent hyper- 
geometric equation, 



d 2 f 



dz 



(B5) 
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which has independent solutions M (£, £, z) and U ((, £, z), 
the confluent hypergeometric functions of the first and second 
kindsi2S In order to avoid the use of complex coefficients (this 
may arise if a single form of the solution is used with both 
directions of the electric field), we give the variable substitu- 
tions and homogeneous solutions Qh for n-type (a<0, 6>0) 
and p-type (a>0, 6<0) cases separately. 

The necessary transformations are obtained by substituting 



fth = (ax + b) f(z) z = 



— (ax + bY 
2a 



(n-type) (B6) 



complicated. The small argument behavior of interest is given 
by 



2^0 sin (7r£) 

1 



1 



r(i + c-0r(0 z«- 1 r(2-or(0 



(Bll) 



The apparent singularity is cancelled by the leading factors of 
(ax+b) in Eqs. (IB6b and ( IB7I . It can be shown after some 
algebra that the functions and derivatives have well behaved 
values at the depletion edge given by 



fi/i = (ax + b) e z f(z) z = 



(ax + by 
2a 



(p-type) (B7) 



into Eq. iB4i . The result of the transformation is Eq. JB5i 
with £=| and 

C = 5-^ (n-type) (B8) 



C = i 



l 



(p-type). 



2aA 2 

In both the n-type and p-type solutions, 

f(z) = 11 M((,^z) + l2 U(t,^z) 



(B9) 



(BIO) 



where 71 and 72 are coefficients determined by applying 
boundary conditions. 

Boundary conditions applied at the depletion edge where 
x=w=—b/a require the evaluation of the special functions as 
the argument vanishes, z— >0. The M (£, £, z) function ap- 
proaches unity as z vanishes. The U (£, £, z) function is more 



fl h (w) = 72 



y/— 2ixa 



1 

2aA' J 



or 



Q h (w) 



f 2na 



-12 



r 1 



1 

2aA 2 



(n-type) 



(p-type) 



(B12) 



(B13) 



and 



dx 



or 



dx 



aji + 72 



= 071 - 72 



— 2ay / 7r 
(-2^) 



2ay / 7T 
^ (2 + 2^A2, 



(n-type) (B14) 



(p-type). (B15) 



The other matching solutions occur for the interface at 
x=0 + and require only the evaluation of the special functions 
for typical arguments with no special considerations. 
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